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Photonic quantum systems are among the most promising architectures for quantum computers. It is well 
known that for dual-rail photons effective non-linearities and near-deterministic non-trivial two-qubit 
gates can be achieved via the measurement process and by introducing ancillary photons. While in principle 
this opens a legitimate path to scalable linear optical quantum computing, the technical requirements are 
still very challenging and thus other optical encodings are being actively investigated. One of the alternatives 
is to use single-rail encoded photons, where entangled states can be deterministically generated. Here we 
prove that even for such systems universal optical quantum computing using only passive optical elements 
such as beam splitters and phase shifters is not possible. This no-go theorem proves that photon bunching 
cannot be passively suppressed even when extra ancilla modes and arbitrary number of photons are used. 
Our result provides useful guidance for the design of optical quantum computers. 



Optical implementations of qubits have played an important role in quantum information science 1,2 . 
Photons exhibit an intrinsic lack of decoherence and are simple to control by standard off-the-shelf 
components. Furthermore, photonic qubits for quantum computation are particularly attractive because 
they can be used to interface to various quantum communication applications 3 . Due to the extremely small 
photon-photon coupling available in existing materials, it was at one point believed that optical qubits could not 
be used for scalable quantum computation. However, it is now understood that the process of photon detection 
itself can lead to effective photon-photon nonlinearities. In particular Knill, Laflamme, and Milburn (KLM) 4 
launched the field of linear optics quantum computing (LOQC) by showing that deterministic single-photon 
sources and high-efficiency single-photon detectors allow the realization of scalable, probabilistic quantum 
computation purely with linear optical elements. This holds in spite of the fact that using linear optics alone, 
the success probability of the nonlinear sign shift gate used in the KLM scheme cannot be improved above 1/4 5,6 . 
Since the original KLM proposal, a number of authors have suggested various simplifications, modifications, and 
optimizations 710 . Although these results formally show that scalable quantum computing is possible, the real- 
ization is very demanding in practice due to the large resource overhead arising from the required non- 
deterministic photon detection events. This becomes particularly apparent when considering beamsplitter-based 
two-photon gates, which are probabilistic in nature. In fact all optical two-qubit gates' 1 , including the promising 
nondestructive CNOT gate 1213 , fail in the case where more than one photon is emitted into the same optical 
output mode, due to photon bunching. One proposal to actively suppress such gate failure events in the polar- 
ization encoding uses the quantum Zeno effect, by coupling to the output light fields to atomic transitions 14,15 . 

While the previous discussion concerned mostly dual-rail qubits, less is known about single-rail photonic 
qubits. Here we address the natural question of whether linear optical elements alone can be used to establish 
quantum interference such that the photon-bunching effect, and thus the probability of gate failure, can be 
eliminated for the case of single-rail encoding. We show that this cannot be achieved in the sense that one cannot 
simultaneously and deterministically implement a linear optical two-photon entangling gate and decouple the 
double-occupancy states (in this work we are not concerned with single-photon nonlocality 16 ). We thus prove a 
no-go theorem for photon bunching suppression via all- unitary linear optics. Our result is complementary to the 
recent proof that the single-photon fraction in any of the single-mode states resulting from purely linear optical 
processing (even including conditioning on results of detections) cannot be made to exceed the efficiency of the 
best available photon source 17 . We stress that our proof applies in the setting of deterministic LOQC; we do not 
address the possibility of simultaneously optimizing the suppression of photon bunching and the success 
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probability of LOQC gates in the single-rail setting. Nor do we 
address here what may be gained by adding measurements or other 
non-unitary operations 18 . 

Results 

Problem formulation. In the standard circuit model of quantum 
computing any unitary transformation on n qubits can be 
decomposed as a product of gates, each of which acts nontrivially 
on at most two qubits, and is the identity on the other qubits 19 . 
Likewise, in the linear optics model, any unitary transformation on 
M modes can be decomposed into a product of "linear optical 
elements," each of which acts nontrivially on at most two modes, 
and is the identity on the other M - 2 modes 20,21 . As linear optical 
operations we consider just passive elements, in particular phase 
shifters and beam splitters, as they generate all linear optical 
elements [i.e., all 2 X 2 unitaries, the group U(2)] when acting on 
the same pair of modes. Moreover, when they are allowed to act on 
overlapping pairs of modes, the linear optical elements generate all 
M X M unitaries [i.e., the group U(M)] 20 . 

Single photon optical qubits come primarily in two varieties: 
single-rail 22,23 , where a qubit is represented by the absence or pres- 
ence of a single photon of fixed polarization in one optical mode, and 
dual-rail 19 , where each qubit is encoded into the presence of a single 
photon in one or the other of two spatial optical modes. Polariza- 
tion qubits 24 — where each qubit is represented by two orthogonal 
polarization modes — are formally equivalent to dual-rail qubits 11 , 
so we will use both interchangeably (see Figure 1). Most LOQC 
proposals, including the original KLM scheme, use dual-rail qubits, 
but there is considerable interest in the single-rail encoding scheme 
as well, as evidenced by the numerous experiments devoted to pre- 
paration of arbitrary states of single-rail qubits 25 29 . In the dual-rail 
encoding single-qubit operations are straightforward while non- 
trivial two-qubits are challenging and are implemented probabilisti- 
cally 4 . The reverse is true for single-rail encoding: it is easy to generate 
two-qubit entanglement deterministically, while single-qubit gates 
can be implemented probabilistically 30 . In both the single and 
dual-rail encodings two-qubit gate failure due to photon bunching 
is a major challenge standing in the way of scalability. Let us next 
explain this in more detail. 

The subspace of interest in the single-rail encoding is spanned by 
the 6-dimensional Fock state basis A(JB where 

a) single-rail encoding b) dual-rail encoding 




Figure 1 | Encoding of optical qubits. (a) Single-rail encoded qubits are 
represented by the absence or presence of a single photon of fixed 
polarization in one optical mode. This encoding allows for the 
deterministic generation of entangled states, but with the caveat that 
single-qubit gates can only be achieved probabilistically, (b) Dual-rail 
encoded qubits are represented by the presence of a single photon in one or 
the other of two spatial optical modes. Formally, polarization-encoded 
qubits are equivalent to the dual-rail encoding due to the basis comprising 
two orthogonal polarization modes, e.g., by defining 1 1)10} = Iff) and 
0)1 1) = IV) where \H) and I V) corresponds to a horizontally and 
vertically polarized photon state, respectively. 



yi=|i,4, 4, «14}l°}> ( la ) 

B={af,4 2 }|0)/V2. (lb) 

Here |0) is the global vacuum, aj, a„ and n t are the creation, 
annihilation, and number operator for mode i, respectively, satisfy- 
ing the standard bosonic commutation relations. The subspace A is 
spanned by the computational basis of single-photon states in each of 
the two modes, while B is the subspace of "bunched" two photon 
states in either one of the two modes (double-occupation states). The 
latter are the problematic states causing computational errors in the 
single-rail encoding. 

To define a basis for the polarization (i.e., dual rail) encoding 
consider the creation operators a] p , where i e {1, 2} is the spatial 
mode index, p e {H, V} is the polarization mode index, and all 
commute for all values of i and p. The ith qubit basis states are then 
1 0) j = a] H 1 0) and 1 1 ) t = flj v 1 0) , so that the computational and double- 
occupation (error) subspaces are, respectively, 

^ = \a\ H a\ H , a\ H a lv , a\ v a\ H , 4y4y}l°>. ( 2a ) 

& = {a%/y/2, a$/y/2, 4,«Jy}|0), MM} (2b) 

It is well known that linear optics is insufficient for generating 
non-trivial (entangling) two-qubit gates in the dual-rail encoding. 
Indeed, this, together with the inherent weakness of nonlinear photo- 
nic interactions such as the Kerr effect, is the raison d'etre of the KLM 
scheme. The reason is a straightforward group-theoretic counting 
argument. Consider m = 2K optical modes. If dual-rail linear optics 
could generate an entangling gate between one pair of qubits, which 
could then be swapped to all other pairs of qubits using beam split- 
ters, the set of linear optical elements would generate the group 
\J(2 K ). However, this contradicts the fact that the set of linear optical 
elements can only generate U(2K) 20,31 . This means that it is imposs- 
ible to generate an entangling gate in the subspace A', irrespective of 
what happens in B' . Another perspective on this is offered by results 
showing the unfeasibility of a linear optics Bell state analyzer, which 
rule out the implementation of a deterministic CNOT gate using 
dual-rail qubits 32,33 . However, these arguments cannot be applied 
in the single-rail case where a beam splitter suffices to generate 
two-qubit entanglement in the subspace A. Since we are concerned 
in this work with deterministic gates we shall not consider the pos- 
sibility of photon-bunching suppression and concurrent enhance- 
ment of the success probability of two-qubit gates in the dual-rail 
case. Instead, we devote the remainder of this work to the no-go 
theorem for deterministic gates in the single-rail case. 

No-go theorem for two modes. We now proceed to prove the 
following no-go theorem: it is not possible to implement a 
deterministic entangling gate between two single-rail qubits while 
at the same time removing the "bunched" two-photon states using 
only linear optics. In this section we assume that only two photons in 
two modes are used; in the following section we allow for an arbitrary 
number of photons and modes. 

The Hermitian beam-splitter generator is 

X= - ^fl|«i +a\a^j, (3) 

and the Hermitian phase shifter generators are, for each mode, 

n l =a\a l , n 2 = a\a 2 - (4) 

Let us also define 
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Then it is easily checked using the identities in the Methods section 
that the set of Hermitian operators {X, Y, Z) is closed under com- 
mutation and moreover satisfies su(2) commutation relations. 
Therefore the beam-splitter generator X and the phase shifter gen- 
erators (combined as Z) generate U(2). The most general unitary 
evolution operator we can then construct according to the Euler 
decomposition is the following composite gate: 



U[oi,p,y,5,e) = A^+Me^^y"^^) , 



(6) 



where we have allowed a more general linear combination of the two 
phase shifter generators than Z, since they can be independently 
tuned. 

We are interested in the (reducible) representation of the unitary 
operator (6) in the 6-dimensional Fock state basis A\JB. We shall 
now show that if one demands that U does not couple A and B, then 
without further assumptions it is not possible to implement an 
entangling gate between the two single-rail qubits. In the AUB basis 
we have 



U{oc,p,y,5,e)=diag(he ia ,e , l i ,e^ +l!) ,e 2ia ,e 2,li 



1 



A, 



diagflAe'V (y+ VV a 
where diag denotes a diagonal matrix and where 



. (7) 



Ar- 



eas € 
i sin e 



A 2 = 



cos 2e 



; sin e 
cos e 

ie i(«+d) sin£ 

e ^+' 5 ) cose 
1 



is/2 



sin 2e 



sin 2e 



1 

is/2 
1 

-is/2 



sin 2e 



sin 2e 



(8) 



B,= 



e i(u+H+v+S) cos 2e e K«+P+iY 



sin 2e 
is/2 



sin 2e 
is/2 



sin 2e 
is/2 
sin 2e 



p2i(a+y) 



s/2 



Mf+y) 



2i(a + &) 



e 2i(P + S) 



Note that U((X,R,y,5,e) consists of three blocks (1, A 2 , and B 2 ) of 
dimensions 1 (vacuum), 2 (one photon), and 3 (two photons), 
respectively. Therefore repeated applications of the composite gate 
17 with different values for the angles y., (1, y, 8, e will still have the 
same triple block structure. Hence our construction is general in the 
sense that we need only consider a single composite gate U. 

Now, the problem is that the two-photon block includes B. We 
would like to decouple these two states from the other four. 
Clearly, this can be done by setting the matrix elements U45, 
l/4 )6 , U 5A and U 6A [i.e., the (1, 2), (1, 3), (2, 1), (3, 1) elements 
of B 2 ] to zero. As is clear from Eq. (8), a necessary and sufficient 
condition for this is 



sin 2e = 0. 



(9) 



This is achieved whenever e is an integer multiple of n/2. If we 
pick an even multiple 2n of n/2 we have 



l/(a,j8,y,<5,mr)=diag(l, (-l)V 



j(u+P+y+S) #i(a+y) e 2, 



If we pick an odd multiple 2n 

u(a,p,y,d,(2n+l)^j 



where 



A.= 



.l)"i e -(P+y) 



1 of n/2 we have 
1 

A 3 

(_l)"fc««+a) 
0 



_ e i(a + /J + y + 5) 

0 
0 



0 
0 

MP+y) 



0 
0 



(10) 



(11) 



(12a) 



(12b) 



In both cases U has a block- diagonal structure where the two- 
photons per mode states are decoupled from the other four states. 
Considering just the 4-dimensional block acting on span(.4), we 
have for the even case 



\n Ma+y) 



diag^l, 

= diag(l,(-l)V (/i+5) )(g)diag( 



.i)" e K«+y) 



(13) 



i.e., a tensor product of two single-qubit phase gates. Similarly, for the 
odd case we have 



1 0 
0 0 

0 (-l)"fcW+iO 
0 0 



0 



(_l)»fc<(«+«) 

0 
0 

= SWAP diag ie ,{rt + s) 



A^+P+y+s) 
®diag( !,(-! 



(14) 



where SWAP = 



a non-entangling gate. Unfortun- 



10 0 0 
0 0 10 
0 10 0 
0 0 0 1 

ately this means that, having decoupled span(B) , the action of U on 
the A subspace is equivalent to a tensor product of two single-qubit 
phase gates, and therefore is no longer entangling. 

The reason for this no-go result is straightforward: we are forced to 
set e equal to a multiple of n/2 in order to prevent coupling to the 
span(B) subspace. Once we do this we are left with only four inde- 
pendent parameters (a, ji, y, 8), which are the parameters of the two 
phase shifters [see Eq. (6)]. The latter can only generate phase gates. 
To circumvent this result would require the use of some type of 
nonlinear optical element (see Figure 2). 

No-go theorem for M > 2 modes. One might wonder whether the 
inclusion of additional modes and photons allows to circumvent the 
no-go theorem. In this section we show that this is not the case. 

Consider a linear optical system with M modes containing an 
arbitrary number of photons, of which only two modes-supporting 
two qubits-are the target of linear optical quantum computation. 
We therefore refer to these two modes as the "computational 
modes." Clearly the choice of which two modes are identified as 
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a) b) 




trivial non-trivial 
two-qubit operations two-qubit operations 



Figure 2 | Two-qubit operations for single-rail encoded qubits. (a) Only the restricted set of trivial two-qubit gates prohibits photon-bunching. 

(b) When using non-trivial two-qubit gates, such as the control-NOT or control-Phase gate, photon-bunching cannot be avoided and thus leads to gate 

errors or probabilistic success rates. 



Let C 



M 



the computational modes is immaterial, so that without loss of gen- 
erality we use the first two modes to this end. The computational 
basis still comprises only those states having either zero or one 
photons per mode in the first two modes, except that now an arbit- 
rary number of photons is allowed in the ancillary modes 3, . . ., M. 

Iljf. ! ^j— | (a'j ^ ' |0) 1 , where the occupation numbers 

{ rij }*|_ l are arbitrary. Then the computational basis states are the set 

Am = {Cm |"i,W2e{0,l}}. The basis set for the problematic subspace 
of bunched states is now Bm = {Cm | «i ,«2 > 2} . We assume that the 
photonic system is initialized in span^^)- Furthermore, we assume 
that the initial state is a product state between the computational and 
ancilla modes, i.e., there is no initial entanglement between these two 
set of modes. Our proof does not apply otherwise, so that initializa- 
tion errors are outside of the scope of this work and will have to be 
dealt with by other means. 

Now let V be some arbitrary passive linear unitary optical trans- 
formation on all, or some subset, of the M modes. For example, V 
could be a product of phase shifters and beam splitters coupling any 
pair of modes. Its most general form is 



( M \ 

= expl i Y, BijaJaA, 



0ij = 0j,e 



(15) 



an arbitrary element of U(M) 20 . The operator V replaces U [Eq. (6)] 
from our earlier two-mode discussion, and we shall show that even it 
cannot be used to decouple the computational basis states Am from 
the bunched states Bm while simultaneously implementing an arbit- 
rary unitary transformation on the computational modes. V has the 
following unitary representation in terms of the M x M matrix v = 
{Vyj, in the basis of bosonic creation operators (see Methods): 

M 



Va] V^=^2 v tj aj , ie{l,...,M}, 



(16) 



Note that V commutes with the number operator, i.e., its only non- 
vanishing matrix elements are between Fock states with equal occu- 
pation numbers. Note also that if we write V = exp(/Q) where 
Q — ^2u=i^'j a i a j' men O|0) = 0 due to the annihilation operator 
a p so that V|0) = jO), and hence (0|V|0) = 1. 

At first sight it might appear that the inclusion of ancilla modes 
could be useful. For example, one might hope that the bunched state 
|200) (for M = 3) can be transformed into the computational state 
(002), where the two photons have been transferred to the ancilla 
mode. However, this approach fails for passive all-unitary linear 
optics. The reason is that we must rule out the reverse process, i.e., 
we must impose the following necessary "don't cause errors" condition: 

V cannot transfer any photons from the ancillary modes into the 
computational modes in a manner which transforms a valid compu- 
tational basis state into a bunched state. 



For example, V is not allowed to take the computational basis state 

fl|J |0) e Am to (nf^ |0) e Bm (ignoring normalization). Thus, 
we impose the necessary conditions for all i S (3,..., M} 



(17a) 
(17b) 



(Otfvia}) |0) = 0, 
(o|fl^y(flj) 2 |o>=o 



Then, 



0 = {0\a\v{a]} \0) = {0\a 2 l Va]v' t VaJv 1; V\0) (18a) 

M 

= Vikvu(o\44 a l v \ o ) ( i8b ) 

k,l=l 

M 

=1^(01^4410)+ v ft v ;! (oi4«!"il°} 

k,l=2 

(18C) 

M 

+ 2^v l7t v,,(0|4ai a I|0) 

k = 2 

= vfi <0|«f «t«i l<>> ( 18d ) 



where to arrive at (18b) we used Eq. (16), and to arrive at (18c) we 
used VjO) = |0) along with the standard bosonic commutation rela- 
tions, and to arrive at (18d) we used the fact that both sums in (18c) 
vanish due to (0|a£ =0. We can thus conclude that v ;l must vanish 
for all i e {3, M}. After we impose the additional necessary 



condition (Oja^f 0 



|0) =0, we find that v j2 must vanish for all 



i e {3, . ..,M}. In addition, it follows from unitarity of vthat then also 
V\i and v 2 i must vanish for all i G {3, M] (see Methods). That is, 
we have proved that the "don't cause errors" condition implies that v 
has a block- diagonal structure: 



v = v c @v a , 



(19) 



where v c is a 2 X 2 block (over the computational modes i,j e {1,2}) 
and v a is an (M - 2) X (M - 2) block (over the ancilla modes i,j e {3, 
. . ., M}). This, in turn, can be interpreted as "don't couple". It states 
that, subject to the "don't cause errors" condition, V cannot couple 
the computational and ancilla modes, and in particular cannot cause 
photons to leak from the ancilla modes into the computational, or 
vice versa. For example, V cannot couple the two computational basis 
states |001) and |100) (M = 3). It is conceptually clear that this 
implies that nothing is gained by the introduction of the ancilla 
modes. However, let us provide a formal proof as well. 
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Assume, as before, that the system is initially in a factorized state 
between the computational and ancilla modes: G^a' 3 , . . . ,a\^j x 
F ( a [ ,(?2 ) 1 0) , where F and G are arbitrary polynomials in the creation 
operators, e.g., F[a\,a 2 ) = Y,m u m*=o c mu 



4)' 



where 



c m,,m2 are arbitrary coefficients. Note that because of the block- 
diagonal structure of v it follows that 

mi,m 2 =0 



mi,m 2 —0 

= F'(a\,a\ 



"'2 

: 1 (20) 



where F' is a new polynomial. Here the key point is that the ancilla 
modes do not appear in F' . A similar calculation reveals that 

V G(aJ, . . . ,4) Vt = G(va\v^ , . . . , Valv^ = G' (a\, . . .,4,), 

where G' is a new polynomial. Therefore the action of V on an 
arbitrary factorized initial state is 



VG [4,..., aUF 4,4 |0) 



(21a) 



VG[a^...,a\ 



V T | |VF(4,4) vt | y l°) (21b) 



c (4,... ,4V f4,4)|o> 



(21c) 

For the computational modes this is equivalent to V F^a\,a^j x 
jO) = (a[,a 2 ^ |0), where the no-go theorem for M = 2 holds. 

Discussion 

We have shown that using linear optics alone it is not possible to 
cancel the photon bunching effect while at the same time implement- 
ing a deterministic universal set of logic gates using single-rail photo- 
nic qubits. In spirit our result agrees with previous theoretical work 
concerning the dual-rail encoding case 4,3233 , showing that passive 
linear optics does not involve particle interactions other than those 
imposed by statistics, and can be understood in terms of classical 
wave mechanics. One approach is to then consider additional non- 
linear operations such as photon-detection or absorption 14 , the Kerr 
effect, or light squeezing, in order to enable universal linear optical 
quantum computing. However, our work leaves open the possibility 
that a more general analysis than we have considered here, in par- 
ticular one which accounts for the possibility of a probabilistic 
enhancement (or even optimization) of both photon-bunching sup- 
pression and linear optical gate fidelity, might circumvent our no-go 
theorem. Nor did we consider measurements of the ancilla modes or 
other non-unitary transformations. For example, one might consider 
a non-unitary transformation with maps the bunched state J 200) (on 
a total of three modes) to the computational basis state 1 002), without 
at the same time mapping (002) back into the computational sub- 
space. In this sense the results presented here are thus a starting point 
for a study of the use of additional resources in linear optics quantum 
computing using single-rail qubits. 

Methods 

Unitary representation of V in the creation operator basis. To prove that the 
representation of V in Eq. (16) is unitary one can use the standard bosonic 
commutation relations along with the Baker-Hausdorff formula 



e *A Be -«A =1 



(22) 



where [ m A, B] := [A, [ m -iA, B]] and [,A, B] = [A, B] , for arbitrary operators A and B, 
and a e C, to establish that v = e~ !t> , where 0 — {(9y} is the orthogonal matrix of angles 



appearing in Eq. (15). Indeed, setting a = i, A— ^2 kj&kj a [ a p B = a l> an d using the 
identity |a|oj,a]j — S^a^ we have: 

M = £M (23a) 



2 A,B] = J20 jk e mi [a]a k ,al] = J2t 

jkm jm 



etc., so that [ m A,B] = Y,f =l [ 0 "Vj- Thus 

00 . m M M 

^-E^E^W-EKW 

m = 0 m -j=l ;=1 

so that Vjj — [e iy ],-,-, whence v = e~ w as claimed. 



(23b) 



(24) 



Proof that "don't cause errors" implies "don't couple". Consider a unitary matrix v 
with the block structure 



x y 
0 z 



(25) 



where x and z are square and y can be rectangular. The unitarity condition v T v = I 
yields 



*t 0 \ fx y 



x l x x'y 



(26) 



so that x*x — I and x t y — 0. By unitarity of v the matrix X cannot be zero, so x~ l = x f 
and hence also xx f = I. Thus xx' ] 'y — y = 0. In the context of our proof in the last 
section, we showed that the "don't cause errors" condition implies that both v,i — v l2 — 
0 for all i e {3, M}. This is the 0 block in Eq. (25). The vanishing of the y block then 
implies that also v u — v 2i — 0 for all i e {3, M}, which is the "don't couple" result. 
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